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r TNTRODUCTION 


The "academic" wind tunnel at NPS is inadequate for most 
tasks, due to the excessive level of turbulence at high 
speed. This tunnel was designed for a potential speed of 200 
Enots in the test section, but this value has not been 
reached yet.This project was undertaken with the purpose of 
evaluating, by a different method from that used for 
designing, the performance of the tunnel in order to 
establish possible causes of the problem. 

The two counter - rotating fans could be the cause of 
the trouble since they were designed using a blade element 
theory, that does not take account of the three dimensional 
effects that are present at medium and high solidities. 
These effects cause changes in the apparent velocity vector, 
therefore the blades were analysed by a " New Blade Element 
Method for Calculating the Performance of High апа 


11 


Intermediate Solidity Axial Flow Fans due to Borst 
[Ref. 1: р.1|. This method determines an induced angle of 
attack which changes the apparent velocity past the airfoil 
and it has been found to be a reliable measure of the three 
dimensional effects.The data required for this method comes 
ош the two dimensional airfoil data, that is available in 
the current literature. This method was used in reference 1 
ІР І0) assuming zero drag, and the predictions showed 
excellent agreement with measurements. The small difference 
between the results can be attributed to the assumption of 
zero drag. The performance of the first and second stage 
blades of the academic wind tunnel is evaluated assuming non 
zero drag, so it is expected that the results will. be 
pealistic. 

A possible cause of the inability to reach the specified 
speed could be excessive losses around the tunnel circuit, 


Euunespally at the first diffeser and the corners. The first 


pu 


diffuser has a high angle of divergence which can provides 
flow separation, further losses and turbulence. The 
evaluation of losses around the tunnel is made by means of 
the method given by RAE & POPE [Ref. 2: p.87], which breaks 
the tunnel down into sections and calculates the losses for 
each one. The four corner losses are evaluated by three 
methods: 


ld considering one third of the т 2. to skin 


friction and two third озона катто ат\ 
empirical relation given b E & POPE. Е 2: EAD. 
was used to account for the different суре си бое 
уапе, 


considering empirical values found by testing and 
3 considering variation of resistance coefficient with 
ds E corner vanes given by PANKHURST 

The first diffuser, since it has variable divergence, 
was broken into four parts and the losses evaluated part by 
part through the method given by RAE & POPE [Ref. 2: p.89], 
this shows that the losses are strongly dependent on the 
diffuser angle. 

After all losses were evaluated , a relation between the 
air inflow velocity to the blades and the power output from 
the blades is found and the non-stall operation envelope is 
determined. The temperature and pressure assumed by LARSON 
in the original blade design [Ref. 4: p.5]; are IES 
р= 2,246 (1bf/ft? ) and, therefore the values of kinematic 
viscos1ity 1 1(p))10t 0 000156 (ft? | вес) апа sound velocity sem 
1160 ft /sec, as presented in Appendix C, will be used for 


evaluation of the Reynolds and Mach numbers. 
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II. EVALUATION OF LOSSES 


EE INTRODUCTION 

As mentioned in chapter one,we will use the procedure 
ОЕ л леа by RAE & POPE [Ref. 2: p.87], that consists of 
breaking the tunnel down into cylindrical sections, 
expanding sections, contracting Sections and corners, іп 
order to evaluate the tunnel losses. 

The sections are numbered in such way that those 
sections with the same characteristics, that is, those that 
use the same loss equation, are kept together. Pais 
facilitates the evaluation by computer. Then, the 
cylindrical Sections are numbered from 1 to 4, the expanding 
ions from 5 to IP) the contracting section and spinner 
take the number 12 and 13 respectively, and finally the 
corners from 14 to 17. 

The Figure 2.1 shows the sections and their designated 


numbers, as they are discriminated below: 


SECTION DISCRIMINATIÓN 
1 test section 
2 fan duct 
3 constant area duct 
4 constant area duct 
5 farst diffuser (sec 1) 
6 й diffuser (sec 2) 
7 М diffuser (sec 3) 
8 ü diffuser (sec 4) 
9 Second diffuser (sec 1) 


10 ü dif rüser (sec 2) 

11 p diffuser (sec 3) 

t2 contraction cone 

15 Spinner 

14 first corner with full size vanes 
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15 осоп сетпет with full size vanes 
16 third corner with full size vanes 
17 fourth corner with half size vanes 

The dimensions of the sections, which are used for 
evaluating the losses, are taken from a WEST COAST RESEARCH 
EO. drawing [Refs. 7,8]. 

Шә, елеп section the loss K will be а drop in static 
pressure divided by the dynamic pressure at that section, as 
given by RAE & POPE [Ref. 2: p.87] and rewritten below 

Ko = K(q/qo), 


K = (pi - p£ ) / q loss at each section, 

Ko = coefficient of loss related to test section, 
pl = inlet section pressure, 

pf = outlet section pressure, 

q = dynamic pressure at section, 

qo = dynamic pressure at test section. 


In terms of diameter we get 


Цой- МЕ (Во Ву where 
Do = test section diameter, 
D = local tunnel diameter, 


mma finally, for energy ratio of the tunnel, 


ERt = jet energy/Summation of circuit losses. 
The losses are evaluated by means of a computer prograin 
written in FORTRAN IV language, which is presented in 
Appendix A, and named LOSS. 
When the evaluation of skin friction is needed, it is 
made by subroutine FRIC in the LOSS program, using the 
equation from RAE & POPE [Ref. 2: p.88], shown immediately 


below, 


IX = 2 logig Rey ) оя (eqn 2217) 


lis 


where A = skin friction coefficient and Rey = Reynolds 
number. 
The Reynolds number at each section 15 related to 
Reynolds number at the test section, as follows: 
Rey = VD/p 


From the continuity equation AV = Ao Vo, and hence 


Rey = (Do/p)Vo(Do/D), from which is found 


Rey = 2.63E 04 Vo(Do/D) . (eqn 2.2) 


B. OGGYLINDRICAB SECTIONS 
Uon Test Section 
The losses for this octagonal section are evaluated 
through the eqn 2.3, (eqn 2.44, Ref. 2, p.88) Тот “ШЕ 


equivalent cylindrical section, 
Kos A(L /De)(Do/De)*, (eqn 2.3) 


where De is the diameter of the circle whose circumferen- 
equals the perimeter of the octagon. 


Then, the dimensions for the test section are: 


Length = 8.00 feet, 
Diameter - 4.36 feet. 
This diameter is the test section diameter (Do) used 
in all the equations related to the test section. 
The losses for the test section at each velocity are 
calculated by the computer program LOSS using the equations 
below, which were obtained from eqn 2.3,from the equation of 


the Reynolds number, and from the dimensions above. 


Section 1: Ko= 1.834A , Rey= 2.63E 04 Мо 
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2. Fan Duct Losses 
This section begins after the first corner with a 
КЕЕ ропа cross section, changes to a cylindrical section 
where the two stages operate and ends in another octagonal 
section. From an original wind tunnel drawing (Ref. 8) we 
have the dimensions for this section and, using the same 
procedure as above (test section), we get the following 
dimensions: 
De = 7.50 ft пее Ft 
The same equations used for test section are used 


here and the equations for the computer program are: 
section 2: Ko- .1448A , Rey= 1.54E 04 Vo 


PRX Comsbanbtexrea Ducts 
ea ses rare "evaluated in the Same manner as for 


Ane" test section . The dimensions for these two ducts are: 


ео, зе ре- 13.37 ft, L= 6.00 ft. 
етой 4. Des 12.73 ft, L-* 6.34 ft. 


The equations for the computer program are: 


ction 3: Коз .0051Х , Rey= .86E O4 Vo. 
meetion 4: Ko= .006/A , Rey= .90E 04 Vo. 


со DIFFUSERS 
ше а јето аге broken down into parts, the first 
(between test section and first corner) into four, and the 
second (between the second and third corners) into three, 
the purpose of which is to make the evaluation of losses 
more accurately. The diffuser divergence angle, which 
Strongly affects the evaluation, is taken for each part as 
the difference between the equivalent small and large 
diameters, divided by the length. 
эе ес топ used comes from RAE & POPE [Ref. 2: p.89], 


and it is 


T7 


Ко= А (A/8tan(a/Z) +. ctantaw = (eqn 2.4) 


where A = (1-(D1/D2)7)(De7 DI)? PPY: 


A = skin friction coefficient for Reynolds 
number given by eqn 2.2, and based 


on the mean diameter, 


и = divergence angle between walls, 
Do = test Section diameter, 

Dl = smaller diameter of diffuser, 
D2 = larger diameter of diffuser. 


1. ¿First Dirffuser 
The diffuser cross section from the test sections 
the first corner varies along its 20 feet of length, and its 
height H and width W are given as a function of the distance 
X (ft), along the flow axis, by the equation 


H = 3.52 + „00440 , (eqn 2.5) 


5.02 + .00036 X? (ft) . (eqn 2.6) 


z 
T 


Furthermore, its cross-sectional area is: 
Area- .828 WH (ft? ) 
The following table gives the local values of H, W, 


Area and De as function of X for this diffuser, where 


De - equivalent local section diameter (ft). 
The evaluation of the diffuser angle is made by the 
following equations: 
-For the divergence angle of the equivalent conical 


Gat fusies 


-For the divergence angle between the walls in the 
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TABLE I 
БӘКЕМЕРЕКЗ ЕШ РЕЕЕТУСЕК 1 


Агеа 
14.63 


15.64 
201159 
22027 
51.94 





K -2arc tan((D2-D1)/2L) . (eqm 277) 
vertical and horizontal planes respectively: 


a, parc Ба к=п) ево) (есп 2.8) 


Pe -2arce tan ((w2-w1)/2L), (eq 2288 
where indexes 1 апа 2 refer to small and large local 
Sections and L is the length between them. 

Шон эст її: ЕФ2 11 15ёг “0 ле equations 2.8 and 2.9 
are used because the angles between the walls are different, 
as shown in Table II 

For a conservative analysis, the larger values are 
chosen in computing the loss coefficient Ko. 

With спе dimensyons from this table and the 
equations 2.4 and 2.2 we got the following equations to use 


with the computer program LOSS 


ща 


TABLE ТИ 
PARAMETERS FOR ЕАСН”ГРАЕТЄ 9 1 171 


Ёс.) DISC EE Ва EE a, (deg. ) a, (deg. ) 


20 4.36 4.46 ОВ 2.2 
әлді 4.46 5.12 3. 60 10728 
2:20 Око 2 босоо ЭТА 20866 
528 ОД В диша ОЈ ЭО ДЕ 





section 35: Kos:530/7 ӨШІР» Rey= 2.60E 04 Vo. 

section o: о оо Веу= 2.40Е 04 Vo. 

secbron 12 Ко-222221 71020222 Rey- 1.97E 04 Мо. 

section 8: оО xD ева, Rey= 1.57E 04 Vo. 
2. Second Diffuser 


For the second diffuser, the values of Dl and D2 are 
substituted into equation 2.7 for evaluation от 3016. 
divergence angle, because the angles between the walls in 


the vertical and horizontal planes have no large difference. 


TABLE NE 
PARAMETERS FOR EACH PART ОЕ ОТЕШ 0 


L Dl D2 


7.40 7.99 
7.88 8.89 
9.89 14.01 





With the dimensions on Table III and equations 2.4 
and 2.2, we got the following equations to use with the 


computer program 0202 
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Beet lene? `: Ko=.14A + .0004 , Rey= 1.50E 04 Vo. 
section 10: Ко=.09Х + .0011 , =Rey= 1.37E 04 Vo. 
Section ll: Ко=.05А + .0030 , Rey= 1.00E 04 Vo. 


D: CONTRACTION CONE 

This section takes the number 12; it is located between 
corner four and the test section. The equation used is from 
RAE & POPE [Ref. 2: p.91], and is 


Ko= „3274 Lc/Do, (eqn 2.10) 
where A = skin friction coefficient for the Reynolds number 
Piven @ bye сам 292 with D=(D1+D2)/2; Lc = length of 


eontraction cone. 


The equations for the computer program are: 
Section 12: Ko=.7338A ‚ Rey= 1.34E 04 Vo. 


КОО SPINNER 
This section takes the number 13, and it is located in 
the fan duct. The dimenSions are shown in fig 2.2 апа the 
loss is evaluated using a method given by Nicolai [Ref. 10: 
Е 11-21). The following equations are for use with the 
computer program: 
Rel = 2.45E-04 Vo, 


A (секет); 

Caon- Cap 000616 | Сат, 

ког 990537 Cdo, 

Сао = drag coefficient, 

Cari friction drag Coefficient, 


where Rel Reynolds Number related to body diameter. 
The dimensions are 

и Length *5-02-0 ft, 

D = Body diameter = 2.0 £t, 


Db 


Base body diameter = 0.554 ft. 


29! 





Figure 2.2 Spinner Dimensions. 


Dx “CORNERS 


The corners are evaluated in three different ways: 
(i) RAE & POPE [ВеЕЁ. 2:. рр. 59, № 
(ii) Bradshaw and Pankhurst [Ref. 5 Ben није 


(iii) Pankhurst and Holder [Ref. 3: p PME 
Ву КАЕ С PORE. tlie) loss equations 


Ко= (0.1% (4. 551 (ово Ва) 2) 052 UE (eqn 2.11) 
where Rn = Reynolds number based on chord vane, 
D = equivalent diameter of inlet свовен ва они 


The dimensions are in Table IV, 
With the dimensions from Table IV and the eqn 2.1 90S 
equations to be used in the computer program are found and 


ae cas follow: 


section l4: Ко= .00690 +.3122/ (19510 Баје 29 
Веу = 1.646Е 03 Vo 

section 15: Ко- .00905“ .4118/(Тов10 Вп)“ “28 
Rey - 1.890E 03 Vo 

section l6: | Кос .00094+ .0428/(logig Кп)2:28 x 
Rey = .609Е 03 \о . | 


| 
sections, Kos .00113* .0515/ (10810 Rn)2: 28 
Реу = —334E 039 05 
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TABLE IV 
PARAMETERS OF CORNERS 


Section Chord (Er) D(ft) 


14 1.042 3252 
T5 1.042 7225 
16 1.042 14.00 
17 UL 13357 





From Bradshaw and Pankhurst, an empirical equation based 
ome tests at Reynolds number 2.0Е+05 апа 1.9E+06, gives the 


moss coefficient as 


Ср = 1.2 (0с/9)7:29 , (бп 2.12) 
where U - local flow velocity, 
с = Vane chord, 


y= kinematic viscosity. 
This equation, when related to the test section velocity 
by means of equation Ко= Ср (Do/D)* yields the following 


equations for use with the computer program. 


section 14: Ко= .0130 Џо“ “25 я 
Section 15: Kos .0160 уо · 5 : 
section 16: | Ко- .0023 Мо "25, 
section 17: Ko= 0032 мо“ "25 


Pankhurst and Holder's method is a graphic one, which 
gives the variation of loss coefficient with the gap/chord 
ratio of the vanes. The tunnel we are dealing with has 
бсшпег - vane sections similar to type b of figure 40 or 
reference 3. Taking the dimensions of the corners and vanes 
iEuu the drawing, evaluating the gap/chord ratio for each 
corner and entering the graph on page 93 of reference 3 


5 


pre value of K ( C is used in Pankhurst 4 Holder) is found, 


До 


and used with the equations below: 


Кк 2 ун/ур 02, (еап 2.13) 
Ко=К (Do/D)*, (eqn 2.14) 
where U = local- velocity at entry of corner: 
УН = drop of pressure across the corner, 
D = equivalent corner diameter at entry, 
K = loss coefficient from the graph. 


Sample calculations 


For corner # 1 we have 
gap / chord s (8.36 - 2.48 )/ 91225 = 0220s 
(ра тр 000007 
K 052955 
Ko = олова 


The losses in all four corners, as given by Pankhurst 


and Holder s method are shonn in p b i em 


Gi. SUMMARY SOR LOSSES 

The losses evaluated by the computer program for 
sections 1 through 17 for a test section velocity of 200 
Knots, are shown in Table VI ; values for velocities from 
100 to 200 Knots are shown in Table VII 

The energy ratio given by Ert is multiplied by r 
order to take into account the losses due to leaks and 
Joints. 

The results presented in Table VII include, for purposes 
of comparison, the losses evaluated by the other two 


methods. At the design airspeed of 200 Knots, Pankhurst 4 
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TABLE V 
CORNER LOSSES BY PANKHURST AND HOLDER 


corner gap/ chord (Do/D)^ K 


L 0.4704 0.0967 
2 0.4256 0.1334 
3 0.4752 0.0094 
4 0.4656 02121113 





Holder's method predicts that 245 HP will be required, while 
Bradshaw & Pankhurst (a later work) indicates a need for 170 
HP. Rae & Pope's method requires 196 HP. The power input to 
the tunnel is 300 x fan efficiency(.85%) = 217 HP. 

The description of the computer program LOSS is 


presented in Appendix D. 
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TABLE VI 
LOSS FOR EACH SECTION OF TUNNEL AT VO-200 KNOTS 


SEC Веу10- 6 A Ko "Loss 


‚88 
TLO 
390 
.04 
S 
2027 
205 
223 
205 
262 
KIT 
255 


test 
CyL 


10298 
0 895 
0.04 
0.05 
2299 

. 42 
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TABLE VII 
ТОП ЛЫ» ВАЭЕНКЕЕР METHODS 


VELOCITY TOTAL LOSSES 
au LOSS (B+ P) LOSS (RIP) LOSS (Pai) 
100.0 22.2610 25.4376 33232752 
120.0 38.0835 43.5826 53.7830 
140.0 60.1842 68.9301 85.2725 
160.0 89.2441 ШОО ЗЕ 126.8912 
180.0 125.2569 143.8930 179.1000 
200.0 170.4057 195.9971 244.5512 
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ІІІ. BLADE PERFORMANCE 


А). INERODUGITON 

The "academic" wind tunnel has two stage counter + 
rotating fans, each fed by an electric motor of 150 hp and 
with a fixed speed of 1200 rpm. The main steps used in the 
design procedure were found in the Larson report [Ref. 4, ] 
which contains the original blade design. The equations and 
the assumptions made during the design indicate that three 
dimensional effects were not accounted for. For this reason 
the method given by Borst, which is based on the blade 
element approach and cascade theory for determining these 
three dimensional effects, will be used to predict the 
performance of the blades. Usefully, this method requires 
data of two dimensional airfoils only, in order to determine 
the force on each section of the blade. 

So, the blade system with all assumptions made during 
the design, (see section B - Larson report), and the Borst 
method, (see section C ), are the material necessary to 
evaluate the performance of the fans. Knowing the lift 
developed at each blade station (eqn. 3.11) one can find the 
torque (eqn. 3.12) at each station. A blade stacion 
defined non-dimensionally by the ratio of the radius of that 
station to the blade tip radius; the root station for 08 
tunnel is at .26/ and the tip statron MASAS Then 
adding the torques at all stations,and multiplying by the 


rotational speed, one gets the power for the blade. 


В. THE EARSON REPORT 

This report describes the procedure that was followed to 
design the blades, together with the assumptions made during 
the calculations. 

The most important assumptions were: 


- 652 fan efficiency for each Базен 
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20 от сечи distribution “ot the” flow 


entering the first stage fan is 


m l2 Cür (rt-r)+Cmt ; (еап 3.1) 


where Cmr = axial velocity at blade root, 
Cmt = axial velocity at blade tip, 
т tip radius, 
= locii adius: 
- Maximum test section velocity equal to 200 
knots. 
- Design lift coefficient equal to 0.64, leading 
to the requirement that Cl- .8 
- Airfoil type МАСА 16 - X12 
This design was based on forced vortex flow, in other 
words, the difference in static pressure across the vortex 
is equal to the difference in tangential velocity (dynamic) 
pressure. This leads us to the loading parameter (solidity x 
lift coefficient) and the mean velocity angle, as one can 


see in Appendix A, whose equations are: 


още ИЗ Z2(A=B)C )/( 1+(Са/С1) В, апа (ean 3.2) 
ВЕ 2122272 ЕВЕ (eqn Эв 
А= Сапј у 4 
В- Сап); В 
ЕО > 
ВЕ ono 


Using these equations and the assumpticns mentiored 
above, the dimensions of the blades were calculated, and are 
given on the drawing of the blades (Ref. 9). The blades were 
divided into 11 (eleven) stations I root to tip. tor each 


station the chord, thickness, leading edge angle, etc... 
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were determined. 


С, THE BORST MEI OD 

This method is based on the blade element approach and 
the vortex theory for determining the three-dimensional 
effects, so that two-dimensional airfoil data can be used 
for determining the resultant force on each blade element. 
By the momentum flow theory, Borst obtained the loading 
parameter (solidity x lift coefficient), as a function 25 
the angles Ву, B5, and f, , that are dependent upon an 
induced angle of attack (a;) at values of constant ite 
angle. 


The loading parameter is: 
Ав 50107. (еап 3.4) 


where i = созр, 

Ы - соза; 

С = Кап ; 

D = БапВо ; 

ЯГ 

Cl-=- lift coefficient, 

Cd = drag coefficient, 

Ві - apparent inlet ane с 

Bo = apparent outlet angle, 

бы angle of the mean velocutv МЕСЕ 
у = асап (Cd/C1), 
g = solidity.: 


Now, based on the equations 


Qi = (81 2 Вэ уу? ; and (eqn 3j. 59 


Bn оог (egm 399 
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іші the two dimensional airfoil data and design parameters 
for the blade, one can iterate the equation 3.4 and get the 
induced angle of attack for a specified test section 
velocity and blade station. 

The angle of attack of the blade at any given station 


Bs. 
СП EOE t Ons (еап 3.7) 


Peete stagger angle (Ref. 9, Ref. 4 , р. 7), which is the 
angle between the chord line of the blade station and the 
line parallel to the axis of rotation of blades.(see fig. 
ЕЕ! ) and 

0- pitch angle (adjustable). 
The Figure 3.1 shows the relation between the variables 


given above. 


D. TORQUE AND POWER EVALUATION 

We will use the Borst method (Ref. 1), to get the 
induced angle for correction of the two dimensional angle of 
attack, and the Larson Report (Ref.4) to obtain all the 
dimensions of the blades and the preliminary design 
Semailtions. 

With this material, we are able to follow the procedure 
outlined below to get to the value of torque and power for 
each fan stage. Mas procedure will ead wus to all 
velocities and angles of each blade, finalizing with the 
value of the torque on each blade section. After that, we 
evaluate the torque for the whole blade, and finally the 
torque and power for the fan Stage. 

PRCCEDURE FOR EVALUATION OF TORQUE AND POWER 

STEP 1) From eqn.3.1 for a given test section velocity 
and n- 1,200 rpm, calculate Wl = ( Cm + UŽ )1/2, 

where U- 2q r n. 
AS one can see, Wl = f( r,Vo) since Ст = f( r,Vo). In 


Appendix A the axial flow velocity (Cm) was developed as 


521. 





Figure 3.1 Angle of Attack Relations. 


шие огу, 

ӨТЕР 2) Calculate by = arc ват (АМА Ron 

STEP 3) From the drawing of the blades (Ref. D 01108 
values of chord are found and then the solidity гои 


blade section is obtained from 

g = B e тапа (eqn 3.8) 
where B = number of blades, c = station chord апа AL 
station radins: 


STEP 4) From the "Aerodynamic Characteristics om m IN 
16 - Series Airfoils" - TN 1546 , (Кеб 12 P UM 
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lift and drag coefficient for each section are available at 
MACH number М = Wl / a , where a = sound velocity at 
temperature of 100 degrees Fahrenheit, and for the design 
lift coefficient of 0.64, which was assumed by Larson (Ref. 
4 ,p.5,6). The angle of attack used to enter the Cl vsa 
ЕС 15 found from eqt.3./, Find CY) Cd and evaluate Уу = 
mime tan(Cd/Cl). 

STEP 5) Assuming an initial value of induced angle, 
evaluate the angles given by equations 3.5 and 3.6. 

STEP 6) Evaluate the right side of equation 3.4 using 
the values found in steps 1,2, and 5; evaluate the left side 
of equation 3.4 using the values found in steps 3 and 4. ТЕ 
both sides are equal, the value of induced angle assumed is 
correct, otherwise, iterate. 

STEP 7) The corrected angle of attack, found by using 
Bae amduced angle, is seen in figure 3.1 and evaluated as 


follows: 
и - а - ај . (eqn 3.9) 


рр) CH corrected angle of attack, recalculate the 
step 4 above and get the corrected lift and drag 
Exsurrcrlent, and evaluate y with these values. 

STEP 9) The mean velocity vector and the lift for each 


blade station are evaluated using the following equations: 


fm) = См /cosB. > (еап 3.10) 
ШІ = 1/2 р ми flc, (eqn 3.11) 
where Li - lift at blade station i per unit length, 
шим Мо Ucoen)oTent at blade station, 
Wn = mean velocity vector at blade station, 


а епа лг Б аде station. 
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STEP 10) By definition, the lift force is perpendi um 
to Wm. For evaluation of the torque it is necessary to get 
the component of the resultant force in the plane of the 
fan, as shown im figure 2-2 1 m 

Ris Гъусоз у, 
Fi = Ет сос (Шон ӨШТІ 


where y = arc жап (Са/С1). МіЕһ Са апа СІ соттеспес е 
evaluated at step 8), апа Li,Ri,Fi are respectively the lift 
force, the resultant force of the lift and drag, апа е 


force in the plane of the fan. 





; Ban Y 


Figure 3-2 Lift Developed at Each Blade Station. 


STEP 11)The force per unit length, Fi, averaged between 
the value at one radial station and next, times the 
difference in radius between these stations, is the force 
used for torque evaluation. This force, times the means 
radius of these two stations is the torque. At this ponte 
is necessary to go through step l to 10 again, Iin orc и 
get Fi for the next station. Then the torque сг 211008 


section between twa consecutive өт 1 011. 
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1 =F 27 + E ) БОГЕ ХЕ, (eqn Br 
where 5 (т - ri) length between stations, 
ШОГ. пи 

and i refers to the ¡th blade station. 

STEP 12) Adding the section torques along the blade, 
multiplying the result by the number of blades and by the 
rotational speed, the power required is found. This power 
provides the test section velocity which was specified at 


the beginning of the procedure. 


Б INTERACTION BETWEEN STAGES 

For evaluating the torque and power for the second 
stage, the procedure given in the previous section can be 
used, but with some adjustments for the different geometric 
EDSracteristics of the blade; for example, the chord, 
stagger angle, etc. . The tangential velocity used is the 
tangential velocity of the. blade at the given station plus 
the tangential velocity of the outflow of the first stage at 
that station. It should be noted that the axial inflow 
velocity is assumed the same as for the first stage. In 
figure 3.4 this interaction can be seen. 

The second stage should be able to accept the tangential 
velocity introduced in the flow by the first stage and add 
the same amount of power as added by the first stage. The 


flow leaving the second stage then has no rotation. 
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Figure 3.4 Е о Velocity Triangles. 
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F.  CONCCUSTON 

Using the computer program BORST and the procedure 
discussed in this chapter, Tables VIII through XX were 
obtained. For all these tables, a pitch angle of 2.4 
degrees was used. 

Ihe Table VIII shows the values of the velocities from 
root to tip, for Vo= 200 Knots, of the first fan Ra 
These velocities are defined as follows: 

О = tangential blade velocity (ЕР M 
Gm = axial flow velocity (tt sec), 


Wl - apparent flow velocity (ft/sec). 


The chord at each station is given in feet. 

The Table IX shows the Reynolds and Mach numbers at each 
station, for the first fan stage and at VO= 200 ҚКпосв ЖИ 
is clear that compressibility effects appear at the outer 
20% of the blade only. The maximum Mach number was less than 
.45 and the values of Cl vs Alpha and Cl vs Cd were obtained 
by linear interpolation between the graphs for M=.3 and 
M=.45. 

The Table X shows the angles between the axial flow 
velocity vector and the velocities relative to the blade, 
the stagger angle, the induced angle, the angle of attack, 
and the lift and drag coefficients. These values are for the 
first fan stage, at each station and at Vo= 200 Knots aie 
values of lift coefficient, are distributed along WEN 
stations of the blade from .662 to .896, around the unit mnm 
value of .8, assumed by Larson. These variables are defined 
in Figures 3.1, 3.4, and the angles are given in degrees. As 
expected, at the root, the induced angle is the largest, 
because the greatest interference between the blades occurs 
there. 

The Table XI shows the angle of attack, lift and drag 
coefficients, all with correction for the induced тсе 
Also shown are the lift, the resultant force and the torque. 


This Table refers to each station of the first Садови 
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uu vos 200 Knots. The values of corrected Cl.(CLCR) are 
ши buwted along the bfade from .358 to ./09, differently 
frome that assumed by Larson, which was a uniform .8 
The variables in Table XI are defined as follow: 
ALPCR = corrected angle of attack (degrees), 
CECR = corrected lift coefficient, 
CDCR = corrected drag coefficient, 
и ЕЕ кстсебрег unit length (IDf/ft), 
DELF = tangential force per unit length (1bf/ft), 
DTOR = torque at middle of two stations (lbf-ft). 

The Tables XII, XIII and XIV, apply to the second fan 
Stage. The variables have the same definitions as for the 
first stage. The Table XII shows the inflow velocity 
relative to the blade (W3) larger than Wl. This is due to 
the tangential flow velocity, which has a finite value after 
the first fan stage. At the root the difference between W3 
and Wl is larger than at the tip, this means that the 
tangential flow velocity at the root is larger than at the 
tip. The Table XIII shows the Mach number less than .45, and 
the interpolation described above was again used. The Table 
Pee snows the Cl varying around .8, that is, from .875 to 
.710 аз for the first fan stage. 

The Table XV refers to the second stage, with the same 
variables as Table XI plus the variable BETA1-BETA4. This 
Table shows the difference between BETA1 and BETA4, which is 
around zero. This means that the flow after the second fan 
stage has no tangential velocity. 

The Tables XVI and XVII, are for the uncorrected values 
of the angle of attack, for first stage at Vo-200 Knots. The 
Table XVI shows the values of Cl around .8 . The Table XVII 
shows values of torque (DTOR) larger than those with the 
corrected angle of attack (Table XI). 

The Tables XVIII, XIX and XX are for the second fan 
Stage without correction of the angle of attack. A 


comparison of the velocities W3 from Table XVIII, and W3 


39 


from Table XII shows that the difference, due во Е 
correction of the angle of attack has affected the outflow 
of first fan stage. The Table XIX shows the variation of 
BETA3 as a result of the variation of W3, just mentioned 
above. The variation of Cl was not so large as to affect 
the value of Cd. А comparison of Tables XV and XX shows 
that the interference between the blades reduces the torque 
by about 1/3 at the root and very little at the tip. This is 
further discussed below. 

In order to visualize the effects of the correction aoi 
the angle of attack, two graphs, each with four curves, were 
plotted; these are: 


I Eggurcj 5 ae Ponce ЦЕ force on the blade is shown 
for each station, at a specified test section үе ба 
of 200 Knots and at a pitch angle of 2.4 degrees. The 
curves are referred to the first and second stages of 
the blades, with and without correction of the angle of 
attack. This figure shows again the effects of the 
interference. The pane See forces have the largest 
difference at the root, and are practically equal at 
the tip. The data were taken from Tables ХТ , XV , 
ХУІІ апа ХХ. 


2 Figure’ 2.61 the 128 of attack is shown for ese 
station of the blade, at а specified test sectam 
velocity of 200 Knots and at a pitch апр ено 
degrees. The corrected angle of attack at the 000 
became negative and, at the tip, that angle is 
practically equal to the uncorrected one. This again 
shows the strong effects of interference at the blade 
root. The curves are referred as above, and the data 
were taken from Tables X , XI F r P T P MN 
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еле ее 


267 


TABLE XVII 


DELL 


38 


39 
37 
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.400 
328. 
узе 


362 
314 


„О6 5 
.845 
222 
43. 
49. 


Bod 
683 
020 


.050 
62. 
64. 


698 
29] 
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Е 
FIRST STAGE, 


ATTACK, 

S 
DELF DTOR 

327260 
4.133 

30.458 
62722 

28.37838 
7.101 

25.482 
Ре 

22.419 
2 2057 

21:0 
152025 

19.364 
17241570 

122122 
18:50:11 

Iv SZ 
20.424 

16.647 
20252023 

14.754 


TAREE ATI 


МЕШОК ЕУ АГ ШОН STATION 
PETHOUISCORNECTION OF ANGLE OF ATTACK, SECOND STAGE, 
VO=200 KNOTS 


STATION CHORD U CM W3 
0.267 Өс 01227700 2321389 
0.300 ОЗИ зола сло 22006 
0.350 Өс ел тоа “12207907 250.296 
0.400 В И 93 475 26389905 
0.450 ООЗЕ 10042 249 0/3 
0.500 © Оо 620 36 559 2955493 
0.600 ROSEO SONA №29 743 331 213 
0.700 ss 52005558 122 377 309040 
0.800 (К АМИ 552 Е 011 1410. 238 
0.900 0.411 424.116 109.145 452.155 
1.000 Ө 1 071207 10272790 (451121 
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O и роса ДИ сети ди СЪРИ Блоа AO UA 
ИВО а СОД оро СЕ о E A EA 
600:0 8/9"0 €6€°%  068:2/ 966 2/ 191721 61/6:6/ 00820 
600'0 665970 690272 066769 020769 1%%89 686:0/ 00120 
60020 7590 81672 O€€°S9 8L7°79 675'69 866799 00970 
60070 06970 %//72 0659:66 27986 858796 %99709 00570 
6000 1490  025'2 060796 999796 <66705 088766 05%70 
6000 9690 458'2 0777165 991706 206'17 790776 00%70 
21070 %/770  09€'€ 0/9797 768797 468707 090775 06670 
61070 €€8°O .Z40'$ 0/6707 05986 066026 756706 00670 
С10:0 25870 79/7 O8S°SE %9067С6 6857902 256787 79070 
42 12 уналу 29415 (ИУІЯЯ 9У144 — €VL3H  NOILVIS 
978'901 -4ЯМ04 05922 = ст атому ношта ана т 
SLONA 00Z=0A ‘ADVIS dqNOO3S “ХОУ11У 40 ЯЛЭМУ 40 ЦОТТОЗЯМОО 100Н1ЇМ 
NOILVLS НОУД ЈУ "4409 ЭУЧАа ЧМУ 1411 ‘SH IONV 


XIX 41ЧУ4 


ТАВЬЕ ХХ 
WITHOUT CORRECTION OF ANGLE OF ATTACK, VO=200 ENOTS 
STATION PECE DELF DOOR 
22:017 90-579 29-799 
2900 31.541 2959055 
23510 51505 221111 
. 400 21102 20112222 
.450 32.748 1222-2153 
251010) 36.370 19.446 
"00 40.413 И. 220 
2275910 46.350 17. T90 
. 800 53.734 161217 
‚900 61.554 ог 20 
.000 02222 Boc UI 


© 


.345 
5.484 
“629 
mM 
.472 
.450 
. 046 
9203 
5912 
2663 


0 
0 
0 
0 
0 
0 
0 
0 
0 
ju 
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IV. TUNNEL PERFORMANCE 


A. INTRODUCTION 

The tunnel performance depends on losses through the 
tunnel circuit and on the fan efficiency. 

The performance will be evaluated to ascertain how the 
loss around the tunnel circuit is related to the velocity at 
the test section, and if the fans can provide adequate 
energy to generate this test section velocity. The energy 
provided by the blade system depends upon the pitch angle 
which must be limited in order to avoid stall. The ріс и 
angle is found for the first and the second stages and 
checked if the flow leaving the second stage has zero 
rotational velocity as required. 

The tunnel performance will be shown through the 


operational envelope of the blade system. 


D v" ОКК Ово о 

As we have seen in CHAPTER TWO the losses of each 
Section depend on the flow velocity, which implies that they 
depend upon of test section velocity 

Relating the losses for each tunnel section to the teem 
Section velocity we can get the total losses as a function 
of this velocity, and then relate the latter to flow energy 
required. 

Using the program LOSS with test section уејос ми 
varying from 100 to 200 Knots, by increments of 20 Knots, we 
get the figure 4.1 which shows the energy required (1) ЭВР-) 


versus test section velocity (Vo) 


C. OPERATIONAL ENVELOPE 

Following the method given in Chapter 3, we find the 
power required from the fans as a function of test ёс ОУ 
velocity, with a number of pitch angles of the (15 


second stage blades. 
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T 
= 
pri 
D 
— 
-0 4 
= у 
j 
„200 ‘pr 
e+e 
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¿OO 
Бо 
i бо 100 /2о ¡Yo 160 180 200 
Vo ( Kuer s] 
Figure 4.1 Energy Required to Overcome Tunnel Losses. 


Since we know the losses around the tunnel in relation 
Test section velocity, we can adjust the pitch angle of 
blades until we get the energy coming into the flow to match 
the losses at same test section velocity. The program BORST 


Was used for this. 


In order to test the sensitivity of the blade flow to 
the approaching velocity profile, we compare the performance 
of the blades in uniform and skewed approach flow. 

Using the axial flow velocity profile, originally used 
by Larson and given by equation 4.1 with the program BORST, 
we got the data shown in Table XXI . This Table presents 
the pitch angle for a specified energy and test section 
velocity, the maximum and minimum pitch angle at that test 
section velocity and its related energy. These values were 
plotted and are shown in Figure 4.2 . This Figure has three 
limit curves, root stall, tip stall and fan efficiency, Gee 
root stall curve limit is the locus of points at үптеп а 
specified test section velocity and pitch angle, the blade 
stalls at the root. The tip stall curve limit is for "Siam 
at the tip. The fan efficiency curve limit is the Тай 
where the operation of the fans is limited by the input 
power and the (assumed) fan efficiency, that is, 300 hp x 
.852 217 hp. This means that the power provided by the 
fan system can not exceed 21/ hp approximately. The pitch 
angles for both stages are equal at the design speed of 
Vo=200 Knots. 


Using the uniform flow velocity given by equation 4.2 , 


we found the data shown in Table XXII . This Table shows 
that the fans can operate without stall, only at a зесрисши 
test velocity from 160 to 1809 ЕВ. These values were 
plotted in Figure 4.3 . This Figure has only tche cuc он 


root and tip stall. The fans stall before reaching the 
maximum power allowed by the fan efficiency. The blades of 
second fan stage are practically stalled at Tenc Troci 

The axial velocity distribution specified for Table 
XXI and Table XXII are developed in Appendix A and shown 


below: 
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Cm шог + 203 Vo , and (едп 4.1) 


Ст 


P369 Vo; (eqn 4.2) 


D. FINAL CONCLUSION 

The large subsonic wind tunnel at NPS, powered by two 
Counter-rotating fans, never achieved the design 
specifications. The reason for the poor performance of the 
tunnel was unknown but believed to be due to either poorly 
designed fan blades or to separation in the diffuser. Being 
easier to analyse, the fan blades were chosen for initial 
study. However, to do this analysis, it was necessary to 
compute the losses in the whole tunnel circuit. 

The results can be summarized as follows: 

1) Most of the losses, about 607 of the total, was in the 
Rest diffuser. The exit end of this diffuser has a high 
value of divergence angle in the vertical plane, which can 
be the cause of a possible flow separation. 

2) When the inflow velocity profile to the fans is assumed 
uniform, the no-stall operational envelope is substantially 
reduced, so there is a much greater likelihood of blade 
Stall. 

3) The variation of the pitch angle affects the 
performance of the fans. The pitch angles for both stages 
should be practically the same in order to have no flow 
rotation after the outlet of the second stage, provided the 
inflow velocity profile is the same as that assumed in the 
design. 

[oan conclusion, although this new blade element method 
used to analyse the fan flow predicts some deviation from 
the original isolated blade analysis, it was found that the 
original design was adequate and, therefore, the tunnel 


problem is most probably due to separation in the diffuser. 
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Полет гоо Envelope for Equation 
4.2. 
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АРРЕМОТХ А 
SPECIFIED FUNCTIONAL FORMS 


CORNER 

Тһе evaluation of the corner loss using the turning vane 
loss function presented in appendix B of reference 6 , 
(p.35) with the flow turning angle equal to 90 degrees, 
leads to the equation given by RAE & POPE [Ref. 2: p.90 ]. 
The equation for the turning vane loss coefficient 1s 

КТУ--1.605670Е-01%1.446753Е-020 -2.570748Е-0467+2. 066207 
Е-06Ф3 -6.335764Е-09%* |, 
МІЕПО- 20 “degrees Fs tound кту= ЕЕ 

The equation for turning vane loss Готова 
ок НО Tn ПРО ери Шек 10810 Кп)2:58), 
where 

Rnr - 500,000 (reference Reynolds number, 
according to ев во poc 
and lo810 5007000: = 2537 
Substituting these values in the equation of the turning 
vane loss function, the equation given by Rae & Pope appears 
as Tol lovs: 
К = 0.15 (2/3 +229 7/10 ON 
K S (0. пое 404657 (los 2D) ши E 
Ко = К (Do/D)^ 

LOADING PARAMETER (Larson) 

The loading parameter used by Larson in his report comes 
from blade element theory, where the effects of drag are 
considered. This is also shown in Osborne [Ref. ll: p.147]; 
by equation 

Vu/Wm- 1/2 С1( 1» Са/С1 соғап (В-а)) : 
where 6 = mean blade velocity angle 
The Larson report (рр. 20120. 
cotan (B a) = tank CEE 
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and rearranging the term Vu/Wm as 
Vu/Wm = (U/Cm - (U-Vm)/Cm) Cm/Wm , and using the axial 
velocity triangles, the following relations are assured: 
U/Cm= tanB,, (U-Vu)/Cm = сапВо ) 
Cm/Wm = COSS 
Substituting these relations into the equation given by 


Osborne, the final equation is 


пати (2 (капва -Бапро )созр „)/(1+Са/С1 сапВ,) and 15 ‘the 
equation used by Larson. 
ООА PLOW VELOCITIES. 


Larson assumed in his report the axial flow velocity as: 


Cite = 2 ie (lle 2001 (21) 
where rt= tip radius, r= station radius, Cmt= axial 
A locity at tip blade, and Cmr = axial velocity at the root 


of the blade. Integrating the product of the cross section 
area of the fan duct at the blades times Cm, from root to 
tip,.and setting it equal to the product of the area times 
the velocity of the test section, we find: 
n ог 65 Сос + 41.037 Cmt. 
ШО equation ot Cm, when r is the radius of the blade 
poot. Cm = Cmr, and 
бое OO (2) 
Пил 1 amd (2), the axial velocity profile for our 
tunnel is found as: 
Еш (02052 ( le ft ) + .303 ) Vo 
UNPEORM AXIAL FLOW VELOCITY. 
This axial profile velocity is used only for purposes of 
comparison. From the continuity equation Ст - Ao Vo / А 


and, taking the known values, we find ЕШ = 903 Vo 
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APPENDIX B 
SAMPLE CALCULATIONS 


Evaluation of the blade station to so 
For first stage 
Following the ‘torque and power' procedure given in 
chapter three, at a specified station the data are as 
follows: 
General data: 
- curves of Cl vsy and Cl vs Cd 20. 
foil МАСА 16-K12 at M=.3 and М=. 3. 
- a=1,160 ft/sec 
- р = .0023 slug/ft? 
y -.000166 ft^/sec 
Б” quu 
сг ns 
- Vo = 200 Knots = 337.55 £t/seec 
- В = 4 blades 
- @ = 2.4 deg. 
АЕ specified station: 
- X = .3 
- с 21 
- В = 40.37 deg. 
SEP) Cm (.2034(1-.3)+.303)х337.55 = 150.34 в 
U = 2x3.1416x.3x3.75x1200/60 mT R E ES 
М1 - (150.342-141.372)1/2 = 206.37 ft/sec 
STEP 2) Бүт = ате Бапа зоос 43.24 deg. 
STEP 3) ос ={(4х1.049)/(2х3.1516;. 57 710100. 
STEP 4) Me 206.37/1, 160-09 185 
a 743.24 - 40.37 + 2.4 = 5.27 deg. 
Entering into curve at M=.3 are found с1=. 090-400 
сат= 015 


y = arc tan (.0157:57 Ма . 
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STEP 


STEP 


STEP 
STER 


STEP 


STEP 


55 assume q. = 5.85 жер. 
Бо MEAM 2: 31.54 deg. 


В 585 - 37.39 deg. 
6) right side -((2xcos?^37.39xcos.988(tan43.24- 
pan iS) cos(3/.39-.988) = 
= о 225 
left siden 8 59365.87 =.5164 
The induced angle assumed is adequate and equal 
Нас сг 
29) 


шэг 0202-0138 пер. 
гэ багт! 

т пап о ОЕ. 31) - 1.397 аев. 
2) wmzaes boa cosel/y . 29 = 189.22 ft/sec 


Li = 1/2x.0023x189.22%x.41x1.049=17.71 1bf/ft 
10) Ri = 17.71 / cosl.397 = 17.71 1bf/ft 
Е2 7 ео (37.99 - 1.397) = 
= 14.33 lbf/ft 


NEXT STATION 


At 


БТЕР 


STEP 
STEP 
STEP 


specified station 
X= .35 
C- -958 ft 


B = 46.67 deg. 
1) Cm = (.2034(1-.35)+.303)x337.55 = 146.90 ft/sec 
АВА Шо 35555 75512000760 -164.93 Ег/зес 
М1 = (146.907+164.932)1/2 = 220.86 ЕЕ /зес 
28 О = ШЕ п пази ОУ О. 90) = 48.31 дер. 
К. 953): (2х3. 14 16х.39х3.75) = .4647 
4) Мэго2209425/ 127 05 „19 


C 
и 


п:48.31 - 46.67 + 2.4 = 4.04 deg. 
Entering into curve at M=.3 are round Cl=.33 and 
Cd = .015 


ое ( 015) 83) = 1.033 агр. 
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ОТЕР Б) assume q;- 4.00 deg. 
Вот 98.31 - 2х4 00 - 40 S Td eg. 


Вз- 48.31 - 4.00 = 44.31 deg. 
STEP 6) right side =((2xcos244.3lxcosl.033(tan48.31- 
-tan40731))/cos(44. 31-1033) == 
= .3864 
left side = .4647/х.83 =. 3866 
The induced angle assumed is adequate and equal 
соу Se: 
SLEP 7) aer 0A А О M 


STEP В) СІ 50555 78200 


Yep атс tan (.01/7455) =e Q 
STEP 9) Wm= 146.9/cos44.3 = 205.25 ft/sec 
Li = 1/2х.0023х205.252х.455х.958-21.11 ДИ cine 
STEP 10) Ri = 21.11 / соз1-.259 = 21.11) 1b не 
Е1-221:117со5(иєг Ї 02259051 
= ad So oul ity) Е 


For: both SCICON =a andae S 
STEP 11) T= (14.33+15.43)/2 x3.75x.05% 225 Ei H: 
ВО ЕВЕ 


After that, adding all station torques and multiplying 

by the rotational speed, the power is found. 
For the second stage, the same procedure сап Бе 
followed with some adjustments as mentioned before in 


chapter) Section E: 
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APPENDIX C 
DESCRIPTION OF COMPUTER PROGRAMS 


The following pages contain the  FORTRAN PROGRAMS 
developed to evaluate the blade performance and the tunnel 
losses. 

The LOSS FORTRAN is a program to evaluate the losses 
around the circuit of the large Academic Wind Tunnel of NPS. 
This program uses one subroutine called SKIN to calculate 
the skin friction when the Reynolds number is given; it does 
this by means of equation 2.1. The main program uses the 
equations for Ko and Rey for each section of the tunnel, 
given in chapter two. The meaning of variables used in the 
program are given at beginning of the main program and 
Subroutine. The listing of the program is presented in 
Appendix D. 

The BORST FORTRAN is a program used to evaluate the 
performance of the blades of both stages. This program uses 
the procedure presented in chapter three section D, to 
evaluate the torque and power of the first stage, and 
section E to evaluate the torque and power of the second ; 
Moses Six Subroutines as follows: FIRST is a subroutine 
to evaluate the parameters of the blade of each station at a 
Specified test section velocity, for the first stage. 
SECOND 15 а Subroutine like FIRST but it evaluates for the 
second stage. CLNEW is a subroutine to compute the lift 
coefficient for each section of the blade using Newton's 
Forward and Backward Interpolation Formula, with the data 
given by curves of Cl уз а in reference 12 .  CDNEW is а 
subroutine to compute the drag coefficient for each section 
using the same formula as in CLNEW, with data given by the 
curve of Cd vs Cl in the same reference as in CLNEW.  INDUC 


is a subroutine used to compute the induced angle of attack 
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through iterations using the equation woe LIFT 15 2 
subroutine to compute the torque and power for second stage 
of blades calling all the subroutines except FIRST; this 
subroutine is called by the main program after it evaluates 
the torque and power of the first stage. 

This program is listed in Appendix E and a description 
of parameters appears at the beginning of either the main 


program or the subroutine. 
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